A high-efficiency quantum non-demolition single photon number resolving detector 
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We discuss a novel approach to the problem of creating a photon number resolving detector using 
the giant Kerr nonlinearities available in electromagnetically induced transparency. Our scheme can 
implement a photon number quantum non-demolition measurement with high efficiency (~99%) 
using less than 1600 atoms embedded in a dielectric waveguide. 
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In recent years we have seen signs of a new techno- 
logical revolution, caused by a paradigm shift to infor- 
mation processing using the laws of quantum physics. 
One natural architecture for realising quantum informa- 
tion processing (QIP) technology would be to use states 
of light as the information processing medium. There 
have been significance developments in all optical QIP 
following the recent discovery by Knill, Lafiamme and 
Milburn that passive linear optics, photo-detectors, and 
single photon sources can be used to create massive re- 
versible nonlinearities^]. Such nonlinearities are an es- 
sential requirement for optical QIP and many commu- 
nication applications. These nonlinearities allow efficient 
gate operations to be performed. In principle, fundamen- 
tal operations such as the nonlinear sign shift and CNOT 
gates have been demonstrated experimentally 0, EI 0- 
However, such operations are relatively inefficient (they 
have a probability of success significantly less than 50%) 
and so are not scalable, due primarily to the current 
state of the art in single photon sources and detectors. 
Good progress is being made on the development of sin- 
gle photon sources 0,0. Absorptive single photon reso- 
lution detection is possible^, IE with efficiencies up 
to ~90% (visible spectrum) and ^30% (infrared, mi- 
crowaves). However, true universal optical QIP will re- 
quire significant further improvements in detector effi- 
ciencies, which will likely require a drastic change of ap- 
proach to detection technology [jl Hoi Hlf . 

In this letter, we propose an implementation of the 
quantum non-demolition (QND) single-photon detec- 
tion scheme originally described by Imoto, Haus and 
Yamamoto[l2], with the required optical nonlinearity 
provided by the giant Kerr effect achievable with AC 
Stark-shifted electromagnetically induced transparency 
(EIT) jljj. We show below that the scheme uses about 
1600 EIT atoms and a weak pulse in the probe mode to 
achieve an error probability less than 1%. The effect of 
the QND measurement in turn means that signal pho- 
tons are not destroyed and can be reused if reauired|l4j. 
Furthermore, for a signal mode in a superposition state 
(such as a weak coherent state) the number- resolving 
QND measurement projects the signal mode into a def- 



inite number state 0], so the detector can be used as a 
heralded source of number eigenstates. We focus on EIT 
as an example here since it has already been used suc- 
cessfully to demonstrate cross-Kerr nonlinearities at low 
light levels in rubidium [lH \v\ . From the perspective of 
realising our detector application, we concentrate on po- 
tential condensed-matter mechanisms of EIT. However, 
clearly any system capable of producing a comparable 
form and strength of Kerr interaction can be used, in- 
cluding optical fibersTSj, silica whispering- gall ery micro- 
resonators^ and cavity QED sy stems QUO]- 

Before we begin our detailed discussion of the EIT 
detection scheme, we first consider the photon number 
QND measurement using a cross-Kerr nonlinearity |l 2l 
l21l . The Kerr Hamiltonian has the canonical form 



H, 



QND 



= hxa ■ ac'c 



(1) 



where the signal (probe) mode has the creation and de- 
struction operators given by a) , a (c^ , c) respectively, and 
X is the strength of the nonlinearity. If the signal field 
contains n a photons and the probe field is in an initial 
coherent state with amplitude a c , the cross-Kerr optical 
nonlinearity causes the combined system to evolve as 
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We observe immediately that the Fock state \n a ) is un- 
affected by the interaction, but the coherent state \a c ) 
picks up a phase shift directly proportional to the num- 
ber of photons n a in the \n a ) state. If we measure 
this phase shift using a homodyne measurement (de- 
picted schematically in Fig. ^ we can infer the number 
of photons in the signal mode a. The homodyne ap- 
paratus allows measurement of the quadrature operator 
x(4>) = ce 1 ^ + c^e - '"^, with an expected result (£(</>)) = 
2Re [a c ] co s <5 + «2Im [a c ] sin 5, where S — </> + n a xt. For 
a real initial a c , a highly efficient homodyne measure- 
ment of the momentum quadrature Y = x(ir/2) would 
yield the signal (Y) = 2a c sin (n a xt) with a variance 
of one, thus giving a signal-to-noise ratio of SNRy = 
2a c sin (n a xt). If the input in mode a is either the Fock 
state |0) or |1), the respective output states of the probe 
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FIG. 1: Schematic diagram of a photon number quan- 
tum non-demolition detector based on a cross-Kerr optical 
nonlinearity [lj. The two inputs are a Fock state \n a ) (with 
n a = 0, 1,..) in the signal mode a and a coherent state with 
real amplitude a c in the probe mode c. The presence of pho- 
tons in mode o causes a phase shift on the coherent state \a c ) 
directly proportional to n a which can be determined with a 
momentum quadrature measurement. 



mode c are the coherent states \a c ) or \a c e txt ). Using the 
momentum quadrature measurement, the probability of 
misidentifying one of these states for one another is then 
terror = ^erfc(SNRy /2 y/2) . A signal-to-noise ratio of 
SNRy = 4.6 would thus give Perror ~ 10~ 2 . To achieve 
the necessary phase shift we require a c sin(xi) ~ 2.3, 
which can be achieved in a number of ways dependent 
upon the range of values available for a c and xt- For ex- 
ample, we could choose a c ^> 2.3 with xt small and sat- 
isfy the above inequality; alternatively we could choose 
Xt = 7r/2 with a c = 2.3. The particular regime chosen 
depends on the strength of the Kerr nonlinearity achiev- 
able in the physical system. 
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FIG. 2: Schematic diagram of a polarization-preserving pho- 
ton number quantum non-demolition detector based on a 
pair of identical cross-Kerr optical nonlinearities. The signal 
mode is a Fock state with an unknown polarization, which 
is resolved into orthogonal polarization states by a polarizing 
beam splitter. The phase shift applied to the probe mode 
is proportional to n a , independent of the polarization of the 
signal mode. 

Figure [21 generalizes the detector shown schematically 
in Fig. n to the case where the polarization of the in- 
put state is resolved into different paths by a polarizing 
beamsplitter. In general we may wish to apply different 
phase shifts to the two distinct polarizations, but in the 
case shown in Fig. [21 an identical phase shift is applied to 
each path, so that the detector is insensitive to the po- 
larization of the input state. This is a particularly useful 
approach when the efficiency of the EIT system and/or 
the optical propagation path (e.g., as provided by a pho- 
tonic crystal waveguide optimized for either TE or TM 
modes) is polarization-dependent. 

We now address the generation of the cross-Kerr non- 



linearity required to perform the QND measurement. 
We consider a model (depicted in Fig. |3J) of the non- 
linear electric dipole interaction between three quantum 




FIG. 3: Schematic diagram of the interaction between a four- 
level M atom and a nearly resonant three-frequency electro- 
magnetic field. We note that the annihilation of a photon of 
frequency uik is represented by the complex number fit,. 

electromagnetic radiation fields with angular frequen- 
cies ui a , ujp, oj c and a corresponding four-level J\f atomic 
system 22] . Mode b should be thought of as a pump or 
coupling field: by choosing the correct conditions, we 
can factor both the coupling field and the atom out of 
the evolution of the atom-field system, creating an ef- 
fective cross-Kerr nonlinear interaction between modes 
a and c. The effective vacuum Rabi frequency for each 
mode is defined as \£lk\ 2 — i&k/VkA) A^ Auik/Str, where 
Ofc = 3A 2 /27r is the resonant atomic absorption cross sec- 
tion at wavelength = 27rc/a;».|23j. r\k is the refractive 
index of the waveguide material, A is the effective laser 
mode cross-sectional area, A^ = /fee 2 w 2 /27reom e c 2 for 
a transition with oscillator strength //., and Aa>& is the 
bandwidth of the profile function describing the adia- 
batic interaction of a pulsed laser field with a stationary 
atomUllilli^. 

It is difficult to achieve a substantial vacuum Rabi fre- 
quency using free-space fields |5^. but encapsulating one 
or more atoms in a waveguide (such as a line defect in a 
photonic crystal structure) allows field transversality to 
be maintained at mode cross-sectional areas on the order 
of A » (A/377) 2 . Consider, then, a two-dimensional pho- 
tonic crystal waveguide constructed from diamond thin 
film (77 = 2.4), with nitrogen-vacancy color centers fab- 
ricated in the center of the waveguide channel [30L l3l| . 
The optical transition at 637 nm in NV-diamond has 
an oscillator strength of approximately 0.12, within a 
factor of three of rubidium, which has been used suc- 
cessfully to demonstrate cross-Kerr nonlinearities at low 
light levels [3 [13 • An EIT transmission window of about 
8 MHz has been observed experimentally [3(3, so a pulse 
with Awfe/27r < 5 MHz should propagate through this 
window with negligible loss. The corresponding vacuum 
Rabi frequency is therefore f2 ~ 3.6 MHz. 

We consider a number N of J\f atoms, fixed and sta- 
tionary within a cylinder that is narrow but long com- 
pared to the optical wavelengths, with the three fre- 
quency modes of the system driven by Fock states con- 
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taining n a , rib, and n c photons, respectively. If the du- 
rations of the three pulse envelope functions are long 
compared to the lifetime of atomic level |2), the evo- 
lution of the amplitude where all of the atoms are in 
the ground state |1) is simply given by |1, n a , rib, n c ) — > 
e~ lWt |1, n a , ribjUc). In general, W is complex (see Eq. 
(139) in Ref. |23); other states contribute to the full 
evolution and there is photon absorption and loss in the 
system. However, for our detector we require the proba- 
bility of even single photon loss from mode a to be very 
small, with a real W preserving the norm of |1, n a , rib, n c )- 
For this, we assume that the laser frequencies io a and 
Wb are both precisely tuned to the corresponding atomic 
transition frequencies, so that the EIT Raman resonance 
condition is satisfied. Then, in the weak-signal regime, 
we assume that the decoherence rate 7^ is dominated 
by spontaneous emission from atomic level \k), and that 
\^a\ ^5 72, so that W is given by 



W = 



N \tt a \ 2 |fi c | 2 n a n c 



u c \Q b \ 2 n b + i (74 |^b| 2 n b + 72 |S! C | 2 n c 



(3) 



where v c = w c — W43. In principle, in this regime we must 
have v c \ fib | 2 rib ^ 74|^b| 2 "b+72|^c| 2 ^c to obtain a nearly 
real W and a low residual absorption. As we shall show 
below, in practical cases where \Wt\ <C 1, this constraint 
can be substantially relaxed. For NV-diamond, 7.7 1 = 
2 x 25 ns 12 8}, [29l . and the spin decoherence lifetime is 
0.1 ms |3JJ, |3l|, so for N < 10000, dephasing can be 
neglected and |^ a |/72 ~ 1- 

Under these conditions, the state \l,n a ,rib,n c ) simply 
acquires a phase-shift which is the basis for the emer- 
gence of the approximate cross-Kerr non-linearity!^ ■ In 
general, when the pump and probe fields arc intense co- 
herent states (parameterized by otb and a c , respectively), 
the evolution of the state |1, n a , ocb, (Xc) can be approx- 



imated as \l,n a ,ab,a c e 



AO 



9 and the residual absorption re is defined for the case 

\n b \ 2 = \n c \ 2 by 



N \SIJU 



\a b \ 2 + i [li \a b \ 2 + 72 |a c | 2 ) 



(4) 



This is equivalent to a damped evolution generated by 
the cross-Kerr Hamiltonian of QJ, with 9 = xt- 

What values of 9 — and, therefore, SNR> — are achiev- 
able? To establish an estimate we need to make several 
assumptions about the physical system and its geometry. 
We assume that the interaction region (where the light 
and Af atoms interact) is encapsulated within the pho- 
tonic crystal waveguide described above, and that the 
pulses have weakly super-Gaussian profiles so that the 
bandwidth-interaction time product is Aujkt ~ 37r, giv- 
ing |f2 a | 2 £ ~ 817772/871. Suppose now that the largest 
phase shift that can be applied in practice (without sig- 
nificantly distorting the signal pulse) is # max , and that 



the corresponding value of a c needed to obtain a given 
signal-to-noise ratio is therefore a c = SNRy/2 max (if 
#max -C 1). Using Eq. (@J, we can calculate the mini- 
mum number of atoms and the corresponding minimum 
detuning v c needed to generate this phase shift, by tak- 
ing the real part and solving for v c explicitly. As v c is 
required to be real we find 
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7 2 SNRy / 74 |a fc | 2 
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with 



\tla\ 2 t 



74 \a.b\ 2 + 72 \a c \ 2 
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(5) 



(6) 



N min and v c 



v cmim we find that 



where the angle 9 



When we choose N 
k = 9, so that the state \l,n a ,a.b,a c ) evolves according 
to \l,n a ,ab,a c e- {1+i ^ 6 ). Therefore, when < 1, the 
residual absorption of a signal photon in mode a can also 
be made intrinsically small, even though the detuning is 
not large compared to the absorption linewidth|32|. 

Figure 01 shows the minimum number of NV-diamond 
color centers as a function of the maximum single-photon 
phase angle for three different values of the error prob- 
ability Pcrror- We note that 74 = 72 for the optical 
transitions in NV-diamond, and for convenience we have 
chosen (rib) = 10(ti c ), requiring a minimum detuning 
fcmin/72 = 1.1 in all three cases. Note that the mini- 
mum number of atoms needed to obtain a given phase 
shift decreases as the phase shift increases, because the 
constraint that the signal-to-noise ratio remain constant 
allows the values of \ctb\ and \a c \ to decrease as max in- 
creases. As an example, we choose P e rror = 0.01 and 
max = 0.01 radians, requiring (n c ) = 5.6 x 10 4 to main- 
tain SNRy = 4.6. Therefore, TV w 1600 is sufficient to 
achieve the desired phase shift, resulting in a residual 
absorption of less than 1%. 

There is considerable flexibility in the engineering de- 
sign parameter space for this implementation of the QND 
detector. For example, if we choose # max = 0.1, follow- 
ing the design procedure outlined above for P 0rrO r = 0.01 
leads to a residual absorption of almost 10% for N w 160 
and (n c ) « 560. However, if we increase the number of 
atoms to 800 and the detuning to v c = II72, then the 
absorption is reduced to 1%. Note also that the detector 
can perform a QND measurement on a Fock state with 
n a > 1 with single-photon resolution. For example, as n a 
increases from 1 to 2, the phase shift n a 9 doubles, and 
the SNR also increases from 4.6 to 9.2 for constant a c . 
The detector sensitivity improves until the phase shift 
becomes so large that one of two fundamental limits is 
reached: either the SNR decreases below the 1% error 
threshold, or the strong nonlinear interaction begins to 
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FIG. 4: Plot of the minimum number of NV-diamond color 
centers needed to generate the phase shift # max for three dif- 
ferent values of the error probability Pcrror in the case n a — 1 . 
We have chosen (rib) = 10{n c ), requiring a minimum detuning 
^cmin/72 = l.X. 



significantly distort the pulse profile of the signal Fock 
state. 

In summary, we have presented a scheme for a highly 
efficient photon number quantum non-demolition detec- 
tor (with single-photon resolution) based on the cross- 
Kerr nonlinearity produced by an EIT condensed matter 
system with approximately 1600 color centers. We have 
explored several different operating regimes, and we have 
examined in detail the performance of the detector for an 
NV-diamond photonic crystal waveguide system. In par- 
ticular, we have shown that efficient detection is possible 
with small phase shifts, which will likely be necessary 
to ensure that the EIT optical nonlinearity doesn't dis- 
tort the pulse envelope of the signal state. Future mod- 
elling will address detector performance for pulsed Fock 
state profile functions, and much experimental work re- 
mains to be done to implement such a detector. For 
example, fabricating EIT atomic or molecular systems 
into a dielectric waveguide is challenging but feasible. A 
method for orienting the color center spins uniformly in 
such a condensed-matter system must be found, and spa- 
tial hole-burning techniques will be needed to overcome 
the effects of inhomogeneous broadening on the trans- 
parency of a condensed matter medium |33j|. The overall 
efficiency of the detector is likely to be limited by the 
efficiency of the homodyne measurement of the phase 
shift, which will depend on the degree to which the ho- 
modyne detector can be spatiotemporally mode-matched 
to a single-photon signal. Nevertheless, EIT provides us 
with the best known candidate mechanism for the imple- 
mentation of the original QND proposal by Imoto, Haus, 
and Yamamoto 12J, and even a weak nonlinearity could 
allow efficient reuse of resources in linear optics quantum 
computation schemes. 
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